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ABSTRACT 

CN ■ 

^ _ Aims. This paper examines the outflows associated with the interaction of a stellar magnetosphere with an accretion disk. In particular, 

we investigate the magnetospheric ejections (MEs) due to the expansion and reconnection of the field lines connecting the star with 
^ ■ the disk. Our aim is to study the dynamical properties of the outflows and evaluate their impact on the angular momentum evolution 

of young protostars. 

Methods. Our models are based on axisymmetric time-dependent magneto-hydrodynamic simulations of the interaction of the dipo- 
f\j . lar magnetosphere of a rotating protostar with a viscous and resistive disk, using alpha prescriptions for the transport coefficients. 

Our simulations are designed in order to model: the accretion process and the formation of accretion funnels; the periodic infla- 

tion/reconnection of the magnetosphere and the associated MEs; the stellar wind. 
(*V^ Results. Similarly to a magnetic slingshot, MEs can be powered by the rotation of both the disk and the star so that they can efficiently 

rjy remove angular momentum from both. Depending on the accretion rate, MEs can extract a relevant fraction of the accretion torque 

and, together with a weak but non-negligible stellar wind torque, can balance the spin-up due to accretion. When the disk truncation 
" approaches the corotation radius, the system enters a "propeller" regime, where the torques exerted by the disk and the MEs can even 

balance the spin-up due to the stellar contraction. 

Conclusions. Magnetospheric ejections can play an imp ortant role in the stellar spin evolution. Their spin-down efficiency can be 
compared to other scenarios, such as the Gho sh & Lambl X-wind or stellar wind models. Nevertheless, for all scenarios, an efficient 
spin-down torque requires a rather strong dipolar component, which has been seldom observed in classical T Tauri stars. A better anal- 
ysis of the torques acting on the protostar must take into account non-axisymmetric and multipolar magnetic components consistent 
1 1 with observations. 
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>: 

"xl" 1. Introduction tracts angular momentum from the star along the field lines con- 

" necting the star with the disk in the region beyond the corotation 

I Classical T Tauri stars (CTTS) are pre-main sequence stars radiuSi where the disk rotates slower than the st ar , Qn t he other 

^ . showing clear signatures of accretion from a suirounding ac- han(i it has been shown that the efficiency of the Ghosh & Lamb 
; cretion disk ^Edwards et alj U99J; | Hartmann et al.j 1998) and me chanism is drastically reduced because of the limited size 

(NJ . election in the form of colhmated jets d Cabrit et alj [1990] of me connec ted magnetosphere (iMatt & 'Pudritzl f2W5^ and 

■ i Burrows etalj | 1996|) . The evolution of their rotation period rep- the dilution of the poloi dal fi eld beyond the co rotation radius 
L; ; resents an interesting puzzle. As soon as they become visi- ( Agapitou & Papaloizou 200a rZanni & Feireiral l2009l hereafter 

. — . ble after the Class 0-1 embedded phases, a relevant fraction p a p er j) 

■ of CTTS appears to rotate well below their break-up lim it, 

U " with rotation periods around 1-10 days (iBouvier et alJl!993h . ° ther solutions are based on the presence of outflows, draw- 
ee ! Besides, their rotation rate appears to be rather constant dur- in S out angular momentum from t he star-dis k syste m, instead 
ing the accreting evolut ionary phases lasting a few million years of transfering it back to the disk. |Shu et al.| dl994j) proposed 
(llrwin & Bouviedl2009h . On the other hand, these protostars are mat an "X-Wind" launched along the open stellar magnetic sur- 
still actively accreting and contracting so that they would be ex- fa ces threading the disk around corotation can extract a sub- 
pected to spin-up at break-up in ~ 10 6 years. Clearly, CTTS stantial amount of angular momentum from the disk before it 
require an efficient spin-down mechanism to explain their ro- is transferred to the star, so as to cancel at least the spin-up 
tational evolution. torque due t o accretion. Even if models of wide-a ngle X-winds 

Since CTTS are kn own to be magnetically active (see m feasible d Anderson et al.| [200l | Cai et al. ||200i, a fully self- 
e.g. Uohns-Krulll 120071 lYang & Johns-Krifll l201lh . different consistent calculation of the disk-outflow dynamical connection 
magneto-hydrodynamic (MHD) mechanisms of angular mo- ' s currentl y missing. 

mentum removal have been proposed. In the iGhosh & Lambl iFerreira et all (120001) investigated a different magnetic con- 
<l 19791) model, originally developed for pulsars, the disk itself ex- figuration, where a magnetic neutral line is formed at the star- 

disk interface due to the cancellation of the stellar dipolar field 

Send offprint requests to: C. Zanni, e-mail: zanni@oato . inaf . it by the disk field. Such a reconnection site has been envisioned to 
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drive massive unsteady ejection events, mainly powered by the 
stellar rotation. These "Reconnection X-winds" provide a very 
efficient spinning down mechanism for early low-mass proto- 
stars (Class and I objects), able to brake a maximally rotating 
initial core down to observed values. However this model has 
not been designed for CTTS and requires a specific magnetic 
topology, not addressed here. 

Stellar winds provide a spin-down torque extracting angular 
momentum along the open magnetospheric field li n es anchored 
onto t he stellar surface (iMatt & Pudritzl l2005bl: ISautv et all 
1201 lh . lMatt & Pudritzl d2008bl) estimated that the wind mass flux 
is likely to be of the order of ~ 10% of the accretion rate in or- 
der to balance at least the torque due to accretion. These stel- 
lar winds would carry the entire ma ss flux typica lly observed in 
T Tauri jets, which seems unlikely (ICabritll2009h . Indeed, such 
a high election effic iency presents a serious energetic problem 
dFerreira et al. I [2006h : since CTTS are slow rotators and their 
centrifugal push is not strong en ough to drive these outfl ows, 
an extra energy input is required. iMatt & Pudritzl ([2005b) pro- 
posed that it could come from the accretion power carried in 
onto the star by the accreting material. But how to transfer this 
power to a sizable fraction of ejected material remains a criti- 
cal issue. It is now cl ear that the required d riving power cannot 
be of thermal origin dMatt & Pudritzll2007l) . On the other hand, 
the push provided by turbulent Alfven waves, such as those ex- 
cited by the impact of the accretion streams onto the stellar sur- 
face, is likely to rem ain insufficient to drive massive stellar winds 
dCranmeiH2 008. 2009). But more importantly, it is quite tricky to 
assume that some accretion energy would be missing (the frac- 
tion that would possibly feed the stellar wind), while still ex- 
plaining the observed UV luminosity. Indeed, it would imply an 
even hig her mass flux onto the star, hence a higher spinning up 
torque dZanni & Fer reira 20 ll]). 

Another class of ejection phenomena is expected to arise 
because of the expansion and subsequent reconnection of the 
closed magnetospheric field lines. The inflation process is the 
result of the star-disk differential rotation and the consequent 
build-up of toroidal magnetic field pressure. This is the same 
phenomenon that bounds the size of the magnetosphere con- 
necting the star with the disk and limits the efficiency of the 
iGhosh & Lambl mechanism. While semi-analyti cal models have 
foreseen the magnetic field e xpansion (see e.g. lAlv & Ku iiners 
1990; Uzdens kv et alj 12002). different numerical experiments 
showed that plasm a ejection is actual l y associated w i th the 
inflation process dHavashi et al .1 119961:" Goodson et al .] [1993 



Miller & S tonelll997l: iRomanova et aflT 2009). Some observable 
properties of this pheno men on have been discussed, for exam- 
ple, by lHartmanB Q009) and Gomez de Castro & von Rekowski 
d201 lb . Besides, iHartmannl d2002l l2009h suggested that this 
mechanism could enhance the angular momentum loss from the 
star-disk system. 

In this paper we present the results of a series of numeri- 
cal MHD time-dependent simulations to analyze in detail the 
energetics and dynamics of these magnetospheric ejections in 
different accretion regimes and evaluate their impact on the an- 
gular momentum balance of the star-disk system. At the same 
time, we are able to include in our models the effects of stellar 
winds. In Section [2] we present the numerical method and pro- 
vide the initial and boundary conditions employed to carry out 
the numerical experiments. In Section[3]we present in detail the 
outcome of a reference case: we first characterize the dynam- 
ical properties of the simulated outflows (Section 13.11 and I3.21 > 
and, subsequently, determine their influence on the angular mo- 
mentum of the disk (Section [33} and the star (Section [3.41 ). In 



Section [4] we study the impact of the disk accretion rate onto 
the dynamics of the outflows and the stellar spin evolution. In 
Section[5]we discuss the outcome of our models making a com- 
parison with other scenarios proposed to solve the stellar spin 
conundrum. In Section[6]we summarize our conclusions. 

2. Numerical setup 

The models presented in this paper are numerical solutions of the 
magneto-hydrodynamic (MHD) system of equations, including 
resistive and viscous effects: 
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This system expresses the conservation of mass, momentum and 
energy and includes the induction equation to describe the evolu- 
tion of the magnetic field. In the system of Eqs. (Q3 p is the mass 
density, u is the flow speed, P is the plasma thermal pressure, B 
is the magnetic field, <J> g = -GM+/R is the gravitational poten- 
tial, J = V x B/4n is the electric current, and rj m is the magnetic 
resistivity, where v m = rj m IAn defines the magnetic diffusivity. 
The total energy density E is defined as 



E = 



y- 



1 



K M B B 

+ p + 

H 2 8tt 



where y — 5/3 is the polytropic index of the plasma. The viscous 
stress tensor t is given by 



t = 77 v 



(Vh) + (v«) t - -<y-u)i 



(2) 



where rj v is the dynamic and v v = rj v /p is the kinematic vis- 
cosity. The anomalous transport coefficients rj m and rj v are as- 
sumed to be of tur bulent origin and parame trized according to 
an a prescription dShakura & Sunvaevlll973l) . The cooling term 

^cool — T]mJ 

/ + Tr(TT T ) /2tj v is included to balance the viscous 
and Ohmic heating, so that the system should evolve adiabat- 
ically, modulo numerical dissipative effects. We employed the 
MHD module provided with the PLUTO codeQ (Mignone et al. 
2007) to solve the system of Eqs. (1). For a prec ise desc ription 
of the employed algorithm, we refer the reader to lPaper j . 

2. 1 . Initial and boundary conditions 

We employ the same initial and boundary conditions, compu- 
tational domain and resolution of the simulations presented in 
Paper I, where a more extensive discussion about the numerical 
setup can be found. For the sake of completeness, we recall here 
the main characteristics of our setup. 

The two-dimensional simulations are carried out in spheri- 
cal coordinates (R, 0) assuming axisymmetry around the rotation 
axis of the star. We indicate the cylindrical radius r = R sin 9 
and the height z = Rco&9 with lower-case letters. We initially 



PLUTO is freely downloadable at http://plutocode.ph.unito.it 
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consider a viscous accretion a-disk surrounded by a rarefied 
corona threaded by the stellar magnetosphere. The Keplerian ac- 
cr etion disk is mod e led af ter the pol ytropic solution presented 
in iKluzniak & Kital d2000l) (see also iRegev & Gitelmanl [2002t 
Umurhan et al. 2006). The density, pressure, toroidal and accre- 
tion speed and kinematic viscosity of the disk are given respec- 
tively by: 



Pd = PdO 
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surfaces anchored inside the Keplerian disk to corotate with it. 
The coronal density is corrected so as to nullify the centrifugal 
acceleration perpendicularly to the magnetic surfaces. 

The computational domain encompasses a spherical sector 
going from the polar axis (8 = 0) to the disk midplane (8 = n/2) 
and from an inner radius R — R+ up to R = 28.6/?*. The domain 
is discretized with a grid of Ng x Nr = 100 x 214 cells. The grid 
is stretched in the radial direction so that the cell sizes satisfy 
the condition AR ~ RA8. Suitable boundary conditions are im- 
posed to satisfy the axial and equatorial symmetries. The bound- 
ary conditions on the stellar surface R = R+ are carefully cho- 
sen to model a perfect conductor rotating with an angular speed 
Q.* so that in the rotating frame of reference the electric field 
E\ a=a — B x (m - £2* x R) — is zero. Besides, the boundary 
is designed in order to absorb the accretion funnels while forcing 
the rarefied plasma of the surrounding corona to have a density 
and enthalpy suitable to drive a light stellar wind. At the outer 
boundary the variables are extrapolated, also ensuring that the 
area of this boundary directly connected to the central star ex- 
erts no artificial torques on the latter. A detailed descripti on and 
discussion of the boundary conditions is given in lPaper 1 



where a v is the anomalous viscosity coefficient, A = 

(5/11 +64/55 a 2 ) , C s = yjP/p is the isothermal sound speed, 

Vk = y/GM+Jr is the Keplerian speed, e — C s /VkI z= o is the disk 
aspect ratio, pdo and Vk* are the density and Keplerian speed on 
the midplane of the disk at R = /?*. The magnetic diffusivity is 
assumed to be proportional to the kinematic viscosity: 



3 v v 
Vm = a m - — , 

2 a v 



(4) 



so that the magnetic Prandtl number is equal to f m = v v /v m = 
2a v /3a m . 

The corona is represented by a polytropic hydrostatic atmo- 
sphere whose density and pressure distributions are given by: 



//M 3/2 2 (JM+ IK+V"- 

where p a o <s pdo is the density of the corona on the spherical 
surface R = R+. 

We model the stellar magnetosphere as a purely dipolar field 
aligned with the stellar rotation axis. Given the flux function *P*: 
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the field components are defined as: 

1 dW± 1 d¥* 
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where and are the stellar radius and the magnetic field 
intensity at the stellar equator respectively. The magnetic flux 
through one stellar hemisphere is equal to: 



<t>+ = 2tiR\ B R (R i ,,8)sm8d8 
Jo 
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2.2. Units and normalization 

We performed the simulations and we are going to present their 
outcome in dimensionless units. We here provide the normaliza- 
tion factors necessary to express the results in physical units, 
taking into account the typical case of a young forming star. 
The stellar radius /?+ is employed as unit length while, given 
the stellar mass and radius, the velocities are expressed in units 
of the Keplerian speed at the stellar surface Vk* = VGM*//?+. 
Taking pd,o as the normalization density, the magnetic field is 

given in units of -JpooV^, time in units of fo = R*/Vk*, ac- 
cretion and ejection rates in units of Mo = Pd0^*Vk*> powers in 
units of £0 = PdO^i^fL and torques in units of Jo = PdoR*V Kir . 
Assuming M* = 0.5M o , =2/? o andp d = 8.5 x 10~ n g cirT 3 , 
we have that: 
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In order to directly provide the characteristic spin-up/spin-down 
timescales, the torques acting onto the star will be divided by 
the stellar angular momentum, expressed in units of 7*o = 
M*/?*Vk*. The inverse of the characteristic braking timescale 
will therefore be expressed in units of 0: 



T* 



= 10 



-6 



PdO 



8.5x10- 



M 



Q.5M e 



-1/2 



\2R Q 



3/2 



yr 



(8) 



The disk surface is determined by the pressure equilibrium 
^d = ^a, while the disk is initially truncated where B 2 /8n = P&. 
In order to minimize initial transient effects due to the differen- 
tial rotation between the disk and the corona we set the magnetic 



2 This normalization differs from the one employed in Paper I, 
where we normalized the stellar angular momentum in units of J t0 = 
PxiR\ Vk*, so that the braking timescale would be given in units of fo. 
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Table 1. Parameters of the simulations: viscosity coefficient a v , 
magnetic resistivity coefficient a m , magnetic Prandtl number 
'Pm - Iv/i/ni: initial viscous accretion rate of the Keplerian disk 
M&. 



10.0 



Simulation 


ar v 


a m 


V 

' m 


M d /M 


CI 


1 


0.1 


6.7 


1.4 x 10~ 2 


C03 


0.3 


0.1 


2 


4.2 x 10~ 3 


C01 


0.1 


0.1 


0.67 


1.4 x 10~ 3 


El 


1 


1 


0.67 


1.4 x 10~ 2 



2.3. The simulations 

Once the initial conditions are normalized, the problem depends 
on six dimensionless parameters: the disk thermal scale height 
e, the equatorial stellar field intensity B^/Bq, the stellar rota- 
tion rate 5* = /?*£!*/ V K *, the coronal density contrast p a o/pdo, 
the viscous and resistive coefficients a v and a m . Except for the 
transport coefficients, the oth er parameters are the same used 
in the simulations of Paper I: e = 0.1, fi* = 5Bq, 5* = 0.1 
and Pao/Pdo = 10~ 2 . Using the standard normalization given in 
Section 12721 this corresponds to a stellar magnetic field = 1 
kG and a period of rotation of the star = 27rfo/£* = 4.65 days 
with a Keplerian corotation radius R co - = 4.64/?*. 

The transport coefficients a v and a m control, respectively, 
the intensity of the viscous torque allowing the disk to accrete 
and the strength of the coupling of the stellar ma gnetic field with 
the dis k m aterial. As discussed, fo r example, in lUzdenskv et alj 
(2002!) and lMatt & Pudritzl d2005al) . the disk magnetic resistivity 
controls the extent of the disk region that is steadily connected to 
the star: since the opening of the magnetosphere is determined 
by the star-disk differential rotation and the consequent buildup 
of toroidal magnetic pressure, a weaker magnetic coupling (i.e. 
a higher a m ) limits the growth of the toroidal field and there- 
fore inc reases the size of the connected region. For example, in 
Paper I we had to assume a value a m — 1 in order to maintain the 
magnetic connection beyond the corotation radius ("extended" 
magnetosphere). 

Since the main aim of this paper is to study the dynamical 
processes associated with the inflation and opening of the mag- 
netospheric field lines, we assume a stronger magnetic coupling 
(a m = 0.1), so that the magnetic configuration opens up closer 
to the stellar surface ("compact" magnetosphere), where these 
phenomena can strongly affect both the disk and the stellar dy- 
namics. Besides, we consider different values of the viscosity 
coefficient a v = 1, 0.3, 0.1 in order to study the evolution of 
star-disk systems characterized by different accretion rates. The 
summary of the cases presented in this paper is given in TableQ] 
W e also in clude the parameters characterizing the reference case 
of Paper I (case El) that will be considered to do some compar- 



3. Star-disk interaction and magnetospheric 
ejections 

Using case C03 as a representative example, we are going to ana- 
lyze the dynamical properties of a interacting star-disk system in 
which the stellar magnetic field is strongly coupled to the accre- 
tion disk. As already pointed out, in such a situation the star-disk 
differential rotation generates a strong toroidal field component 
and, due to its pressure, the magnetic structure relaxes by in- 
flating and opening the initial dipolar configuration close to the 
truncation radius. The overall picture illustrating the outcome of 
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Fig. 1. Global view of the star-disk interacting system. A loga- 
rithmic density map is shown in the background. Poloidal speed 
vectors are represented as blue arrows. The dotted line marks 



the Alfven surface, where u p = B p / -\J4np. Sample field lines are 
plotted with white solid lines. Thick yellow field lines, labeled 
as (a), (b), and (c), delimit the different dynamical consituents of 
the system indicated in the figure. The image has been obtained 
by averaging in time the simulation outcome over 54 stellar pe- 
riods. 



this process is given in Fig. Q] Four groups of field lines can 
be distinguished: (1) the field lines steadily connecting the disk 
with the star below the magnetic surface (b); (2) the open field 
lines anchored on the surface of star at latitudes higher than the 
position of the surface (a); (3) the open field lines attached to 
the accretion disk beyond the surface (c); (4) the field lines en- 
closed between surfaces (a), (b) and (c) connecting the disk with 
the star, periodically evolving through stages of inflation, recon- 
nection and contraction. An example of this periodic process is 
represented in Fig. [2] The periodicity of these phenomena cor- 
responds to about 2 stellar rotation periods. Anyway, since the 
reconnection processes which are involved are driven by numer- 
ical resistivity, this periodicity has to be considered cautiously. 
On the other hand, this almost periodical behavior allows us to 
use time averages to characterize the long term evolution of the 
system and smooth out transient features. For example, Fig. Q] 
has been obtained by averaging snapshots over 54 rotation pe- 
riods of the star: notice that, because of the time average, the 
fleeting reconnection phenomena are not visible. 

Different dynamical processes are associated with the four 
groups of field lines. In the region inside the magnetic surface 
(b), ("connected disk") the star and the disk can directly ex- 
change angular momentum, the disk is truncated and the ac- 
cretion curtains form. Notice that this region extends within the 
Keplerian corotation radius R co = 4.64/?* so that, beyond this 
radius, the disk and the star do not have a direct magnetic con- 
nection. Therefore, the Ghosh & Lamb scenario cannot be di- 
rectly applied: the disk region beyond corotation, which rotates 
slower than the star, cannot exert any direct spin-down torque 
onto the star. Three classes of outflows correspond to the other 
groups of field lines. A stellar wind flows along the open mag- 
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Fig. 2. Temporal evolution of the periodic infiation/reconnection process which characterizes the dynamics of the magnetospheric 
ejections in case C03. We show logarithmic density maps with sample field lines (white solid lines) and poloidal speed vectors (blue 
arrows) superimposed. The yellow solid lines follow the evolution of a single magnetic surface showing clearly the dynamics of the 
phenomenon. Time is given in units of rotation periods of the central star. 



netic surfaces anchored at high stellar latitudes. A disk wind can 
be accelerated along the disk open field lines, but only the field 
lines closer to the star are characterized by a field strong enough 
for this outflowing component to play a relevant role, as it will 
be shown in Section [33] Finally, the inflation at mid-latitudes 
of the dipolar field lines is very dynamic and it is accompanied 
by outflows that can in principle extract mass, energy and an- 
gular momentum both from the disk and the star. On a relatively 
large scale (10-20/?+, see Fig.|2]i, these ejections detach from the 
magnetosphere in a reconnection event and continue their prop- 
agation as magnetic islands disconnected from the central part 
of the star-disk system, in between the open magnetic surfaces 
anchored into the star and those anchored into the disk. In the 
following we will refer to this type of outflow associated with 
the process of infiation/reconnection of the magnetospheric field 
lines as magnetospheric ejections (MEs). 

The different regions outlined in Fig. Q]can be also charac- 
terized by the amount of poloidal magnetic flux which partici- 
pates in each of them. In Table [2] we show for all the discussed 
cases the fraction of magnetic flux O = J B p ■ dS which crosses 
each region relative to the total stellar flux through one hemi- 
sphere, Eq. (O: we display the stellar wind flux <t>sw, also equal 
to the open magnetic flux of the disk wind, the MEs flux <J>me, 
the flux of the connected disk Ocd and the flux contained in the 



Table 2. Fraction of stellar poloidal magnetic fluxes threading 
different regions of the star-disk system in different cases: stel- 
lar wind <I>sw, MEs <E>me, steadily connected disk <I>cd, magne- 
tospheric cavity <f>MC, accretion columns <1>ac- Also shown are 
the fractional surface of the accretion columns S ac, the posi- 
tion of the truncation radius R t , the anchoring radii at the disk 
midplane of the outermost closed (R cm ) and the innermost open 
(R om ) magnetic surfaces. Different estimates have been done for 
the accretion and propeller phases of case C01 (see Section l4~2l ). 





CI 


C03 


C01 (acc.) 


C01 (pro.) 


El 


fsw/f* 


0.052 


0.066 


0.091 


0.109 


0.091 




0.023 


0.035 


0.039 


0.053 


0.009 


Ocd/®* 


0.269 


0.180 


0.135 


0.102 


0.239 


Omc/®* 


0.656 


0.719 


0.735 


0.736 


0.661 


®ac/®* 


0.194 


0.135 


0.087 


0. 


0.152 


S ac/S * 


0.110 


0.075 


0.049 


0. 


0.088 


RJR* 


1.9 


2.7 


3.0 


4.4 


2.5 


R<m/R* 


2.7 


3.7 


3.8 


4.7 


11.9 


Rom/R* 


3.1 


4.7 


4.2 


5.7 


15.1 



magnetic cavity inside the truncation radius Omc- Obviously we 
have Osw + ^me + <I>cd + ^mc = Since not all the field 
lines steadily connecting the star to the disk are mass-loaded to 
form the accretion funnels, we also show the amount of mag- 
netic flux threading the accretion columns <5ac, clearly equal to 
a fraction of <&cd, and the corresponding surface covering frac- 
tion S ac /Si,. To give an indication of the size of the different in- 
teraction regions we also provide the position of the disk trunca- 
tion radius R t , the anchoring radius R cm of the outermost steadily 
connected magnetic surface (labeled as (b) in Fig. [TJ and the an- 
choring radius R om of the innermost open magnetic surface (la- 
beled as (c) in Fig. [TJ at the disk midplane. Notice that, while 
the magnetic flux is frozen into the stellar surface, its distribu- 
tion can change on the disk midplane. Besides providing a clear 
indication about the relative importance of the different dynami- 
cal components, these quantities can be directly compared to the 
predictions of other models, such as the X-Wind. 

3. 1 . Dynamical properties of magnetospheric ejections 

We are going to characterize the properties of the MEs by in- 
specting their mass, angular momentum and energy fluxes. By 
defining a surface 5 perpendicular to the poloidal flow these are 
respectively defined as: 



M 



■X 
-X( 
-XKI 



dS 



rpu^Up - 



-pu 2 + 



rB^Bp 






j-dS 


4n 




yP 


GMi, 


7-1 


R 



Up + E x B\ 
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(9) 



(10) 



(11) 



As already pointed out, this type of ejection can exchange mass, 
angular momentum and energy both with the star and the disk. 
We can therefore define different contributions to the budget of 
the MEs. Taking for example the mass flux, we can estimate the 
stellar mass input as: 



ME.s 



AnRi 



piiR sin 8 d6 



(12) 



where 0. d and 9 are the anchoring angles onto the stellar surface 
of the lowest open and the outermost steadily closed magnetic 
surfaces respectively, labeled as (a) and (b) in Fig. Q] The disk 
contribution can be calculated as: 



M 



ME.d 



= -4tt f ° 
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rdr 



(13) 
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Fig. 3. Temporal evolution of mass fluxes of the different accre- 
tion and ejection phenomena present in the system: mass accre- 
tion rate measured at the stellar surface {solid line), total mass 
outflow rate of magnetospheric ejections {dot-dashed line), mass 
flux fueling the MEs coming from the star only {long-dashed 
line), stellar wind outflow rate {dashed grey line), disk wind 
mass outflow rate {dotted line). The image refers to case C03. 
Time is given in units of rotation periods of the central star. 

where and R out are the anchoring radii at the disk surface 
of the outermost steadily closed magnetic surface and the inner- 
most open field line threading the disk, labeled as (b) and (c) in 
Fig. Q~l whose footpoint radii at the disk midplane are R cm and 
^ om respectively (see Table 0. These radii are marked in the 
panels of Fig. l7land[T0l The + superscript indicates a quantity 
evaluated at the disk surface H (r) = 2C s /£2k, proportional to the 
thermal heightscale; the poloidal vectors labeled with a + are de- 
fined as the components perpendicular to the surface of the disk: 



for example, in the case of the speed, ul 



-u z + u r H' 



H(r) 



uq (sin 8 + cos 6H') — ur (cos 6 — sin 6H r )\^ r y The total mass 
flux of the MEs can be evaluated by choosing a surface crossing 
both magnetic field lines (a) and (c); for example, by selecting 
a spherical zone with a radius R > 6R+, so that it crosses both 
magnetic surfaces (a) and (c), we define: 



f 

J ft. 



A^ME.tot = 4nR pu R sin 0d6, 



(14) 



where 6> a and 6 C are the angles at which magnetic surfaces 
(a) and (c) intercept the sphere with radius R. When calcu- 
lating the time evolution of the previous fluxes, the anchoring 
angles and radii can vary with time. In a steady situation or 
taking into account time-averaged quantities for a system that 
evolves periodically (as in the considered case), we have that 
A^ME.tot = Mme,s + Mme4- Analogous expressions can be de- 
rived for the angular momentum {Jme.s, Jmea, Jme,^) and en- 
ergy fluxes (£ M e, s , E MEtd , £ M E,tot) by integrating Eq. ([TOj and 
( fTTb on the same surface elements employed to define the mass 
fluxes. Notice that we exploited the midplane symmetry of our 
simulations so that the flux integrals are related to two-sided out- 
flows. 
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Fig. 4. Temporal evolution of the specific angular momentum 
carried by different accretion and ejection component of the sys- 
tem in case C03: specific angular transported by the accretion 
flow through R out (upper panel); total specific angular momen- 
tum carried by the MEs (second panel); specific angular mo- 
mentum extracted by the MEs from the disk only (third panel); 
specific angular momentum of the stellar wind (lower panel). 
Dashed lines show the temporal averages of the quantities over 
the plotted lapse of time. Time is given in units of rotation peri- 
ods of the star. 



mass accretion rate measured onto the surface of the star, defined 
as: 



r*/2 

I pu r sin 6 d8 . 



(15) 



After an initial transient lasting around 20 stellar rotation peri- 
ods, the total outflow rate of the MEs, calculated using Eq. (fT4l > 
with R = 1R*, regularly oscillates around a value corresponding 
to a: 18% of the accretion rate at the stellar surface. This corre- 
sponds to * 15% of the disk accretion rate measured at R out . The 
oscillations correspond to the periodic inflation and reconnec- 
tion phenomena. The stellar contribution to the MEs, Eq. (fT2l . 
is of the order of a few percent of the accretion rate. This means 
that MEs are essentially mass-loaded from the disk and their in- 
ertia is dominated by material coming from the accretion disk. 
The stellar and disk wind mass fluxes plotted in Fig. [3] will be 
defined and discussed in Section l3~2l 



3.1.1. Mass fluxes 

In Fig. [3] we plot the temporal evolution of the mass-loss rates 
of the different outflowing components and compare them to the 



3.1.2. Angular momentum fluxes 

Since the MEs are magnetically connected both to the star and 
the disk, they can potentially extract angular momentum from 
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Fig. 5. Radial evolution of the total specific energy carried by the 
MEs in case C03. The total specific energy (solid line) is given 
by the sum of Poyinting-to-mass flux ratio (dot-dashed line), ki- 
netic energy (long-dashed line), specific enthalpy (triple-dotted- 
dashed line) and potential gravitational energy (dotted line, plot- 
ted in absolute value). The specific poloidal kinetic energy is also 
plotted (dashed line). The slow- and fast-magnetosonic surfaces 
are marked by a vertical line. The plot starts from the cusp of the 
innermost magnetic surface that steadily connnects the disk and 
the star (labeled as (b) in Fig. [1). The figure has been obtained 
by time averaging the energy and mass fluxes over 54 rotation 
periods of the star, from time t — 38 up to t — 92. 



both. In the second panel from the top in Fig. [4] we plot the tem- 
poral evolution of the normalized total specific angular momen- 
tum carried by the MEs, defined as: 



7ME,tot - 



■/ME.tot J ME.t( 



(16) 



while in the third panel from the top we have evaluated the spe- 
cific angular momentum extracted by the MEs from the disk 
only: 



7ME,d - 



ME.d 



Mi 



(17) 



ME.d 



Plots are in units of -\JGM+R+. Both quantities regularly oscil- 
lates in time but, if we look at the time averages, we clearly see 
that the specific angular momentum extracted from the disk is 
around ./ME,d ~ 7 yjGM+R*, while the total angular momentum 
of the MEs is approximately ./ME.tot ~ 10.5 a/GM^RT. Clearly, 
the total angular momentum of the MEs is larger than the an- 
gular momentum extracted from the disk only, meaning that a 
substantial fraction comes also from the star. In Section [331 and 
13.41 we will discuss the effects of these torques on the angular 
momentum evolution of the disk and the star respectively. 



3.1.3. Energy fluxes 

Concerning the energy budget of the MEs, the power extracted 
from the disk only corresponds to ZsME.d ~ 0.2GM*M acc . s //?i n * 
0.05GM*M acc >s /7?*. It is important to point out that the mechan- 



ical power released by the material accreting from R oM down to 
R,„, defined as 



-M m 



" 2 „ 



(18) 



is sufficient to power the part of the MEs coming from the disk. 
The power extracted from the star is comparable and largely de- 
termined by the Poynting flux associated with the spin-down 
torque 7me,s- The enthalpy flux, needed in our simulations to 
give the initial drive to any stellar outflow, is a small fraction 
(« 0.01GM*M acc s //?*) of the accretion power. In case C03, the 
MEs operate as a magnetic sling, powered both by the stellar and 
disk rotation. 

In order to have a better understanding of the asymptotic 
properties and the acceleration efficiency of MEs, we can in- 
spect the energy conversion along the flow. In Fig.|5]we plot as a 
function of the radial coordinate R the evolution along the flow 
of the total specific energy of the MEs: 



<?ME,tot 



Mi 



(19) 



ME.tot 



This plot has been obtained by averaging the mass and energy 
fluxes over 54 stellar periods and it starts from the cusp of field 
line (b), located at * 5/?*. In a stationary situation this definition 
corresponds obviously to an average of the Bernoulli invariant 
over a section of the outflow; 



1 



-u + h 



R Anpu p 



(20) 



given by the sum of kinetic energy, enthalpy /; = yP/(y - \)p, 
gravitational and magnetic (Poynting) energy. On the other hand, 
the more general definition Eq. ( TT~8b allow us to define an energy 
conversion efficiency in the case of non-stationary MEs. We can 
clearly see that the total energy is not conserved along the flow: 
this is mainly due to to the dissipation of the magnetic energy 
due to the reconnection events. Even if the reconnection is con- 
trolled by numerical dissipation and therefore is not physical, it 
leads to a temperature increase, visible at lower radii in Fig. [5] 
The outflow cools down subsequently: since we did not include 
any realistic cooling function, we simply limited the maximum 
specific entropy (P/p y ) that can be attained by the outflow in 
order to preserve code stability. Notice that the kinetic energy 
stops increasing after ~ 10^*: this reflects the fact that, after 
the plasmoids have detached from the inner magnetosphere, the 
acceleration process stops and the propagation of the outflow be- 
comes ballistic. These aforementioned effects limit the accelera- 
tion efficiency of the outflow, whose terminal speed is around 0.5 
Vk*. Finally, consistently with their ballistic propagation, MEs 
cannot be self-collimated thanks to magnetic stresses, as usu- 
ally envisaged for magnetically driven outflows. This confers the 
characteristic "conical" s hape to these magnetosp heric outflows, 
as already pointed out by Roman ova et alJ (120091) . On the other 
hand, they can be in principle confined by some external agent, 
as it will be discussed in AppendixlAl 

3.1.4. Forces 

To complete the analysis of the dynamical properties of the MEs 
for the reference case C03, we take into account the forces that 
drive these outflows. Since mass is accelerated both from the 
disk and the star, we consider an inflating field line anchored in 
the disk at R = 4.3.R* still connecting the disk and the star: we 
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Fig. 6. Forces acting along a field line connecting the star and the disk that is mass-loaded by the MEs in case C03. The left panel 
refers to the forces acting from the disk midplane up to the cusp of the magnetic surface, while the right panel shows the forces 
from the stellar surface up to the cusp. We selected a magnetic surface anchored at 4.37?* at the disk midplane. The plots have 
been obtained by time averaging the forces over 54 stellar rotation periods. We plot the thermal pressure gradient (solid line), the 
centrifugal term (dot-dashed line), the Lorentz force (dashed line) and the gravitational pull (dotted line). Vertical lines mark the 
position of the disk surface (DS) and the slow-magnetosonic point (SM). 



plot the component of the poloidal forces parallel to this field 
line from the disk midplane up to the field line cusp (left panel 
in Fig. [6]) and from the stellar surface up to the cusp (right panel 
in Fig. |6]l to analyze the driving mechanism of the MEs from 
the disk and from the star respectively. The forces have been 
obtained from a time-averaged snapshot, in order to smooth out 
transient features. 

In the left panel of Fig. [6] we show that the acceleration of 
the mass of the MEs coming from the disk is largely due to a 
combination of centrifugal and magnetic effects, as in a typical 
disk-driven outflow. In addition, the thermal pressure gradient 
-V||P is comparable to the centrifugal acceleration at the disk 
surface (DS, defined as the point where the Lorentz force J x B\\ 
changes sign) and crucially contributes to the outflow acceler- 
ation. As it will be shown in Section 13.31 this enhanced ther- 
mal pressure gradient is due to the push of the accretion flow 
against the magnetospheric wall and is responsible for the high 
mass load of the MEs coming from the disk. Thermal effects are 
important since the centrifugal term pu^jr\\ is not sufficient to 
counteract the gravitational pull pg\\ at the base of the flow: as it 
will be more extensively discussed in Section [3~3l in the disk ac- 
celeration region of the MES, the disk rotation becomes strongly 
sub-Keplerian. 

The right panel of Fig. [6] shows that the pressure gradient 
provides the initial thrust of the mass of the MEs coming from 
the star. Even if this term is most probably not of thermal type, 
as it is assumed in our simulations, it is needed to drive any kind 
of stellar outflow from slowly rotating stars, where magneto- 
centrifugal effects are not sufficient to give the initial push. The 
sudden change in the profile of the forces at c s 1R+ happens 
when the material accelerated from the star comes across the 
mass coming from the disk. The latter is characterized by a 
higher density (see, for example, the change in the centrifugal 
push and gravitational pull), confirming the fact that the MEs 
inertia is dominated by the mass loaded from the disk, as al- 
ready discussed in Section [3.1.1l The profile of the Lorentz force 
has some interesting features. We first recall that the Lorentz 



force parallel to a field line in the poloidal plane is related to the 
toroidal component of the force, according to the relation: 

(/ x B) Bp = - (/ x B) Bt , 

which shows that a Lorentz force accelerating (braking) along 
a field line also accelerates (brakes) in the toroidal direction. 
Therefore we can see that the mass loaded from the star is sub- 
ject to a toroidal acceleration close to the stellar surface while, 
getting closer to the part of the MEs coming from the disk, be- 
tween 5R+ < r < 1R+, it is spun-down. This clearly indicates 
that the star is trying to spin-up the plasma attached to this field 
line, therefore loosing angular momentum, countering the ma- 
terial of the MEs coming from the disk that is trying to spin it 
down. 

3.2. Dynamical properties of stellar winds and disk winds 

In our simulations, a stellar wind is accelerated along the open 
field lines anchored into the stellar surface. Due to the energetic 
limitations illustrated in the Introduction, we limit the mass out- 
flow rate of the stellar wind to a few percent of the mass accre- 
tion rate, so that the energy needed to drive initially the outflow 
corresponds to a small fraction (less than 10%) of the power dis- 
sipated by accretion onto the stellar surface. In Fig. [3] it is also 
plotted the temporal evolution of the stellar wind mass outflow 
rate, calculated as: 

M sw =4t:RI pu R sin 6d6 , (21) 
Jo 

where 8 a corresponds to the anchoring angle of the last open stel- 
lar magnetic surface: the outflow rate corresponds on average, to 
w 1.6% of the mass accretion rate measured onto the surface of 
the star (Eq.[T5ll. Comparing the mass outflow rate of the refer- 
ence case with the stellar wind of case El, we find that the mass 
loss rate of case C03 corresponds on average to ~ 30% of the 
outflow rate of case El, despite in the two cases we assumed ap- 
proximately the same injection speed, density and temperature. 
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The mass ejection rate difference is therefore due to a different 
size of the launching area: Table [2] shows that in cases C03 and 
El the closed magnetosphere contains approximately the same 
amount of stellar flux (<E>mc + ^cd), even if in case C03 this flux 
is compressed in a much smaller region closer to the star (see 
R cm in Table [2}; the presence of the MEs, which are almost ab- 
sent in case El, reduces the amount of magnetic flux and stellar 
surface which is available to launch the stellar wind. This ex- 
amples shows that a self-consistent model of stellar winds from 
accreting protostars must take into account the presence of the 
accretion funnels, which can strongly affect the geometry of the 
launching region and the morphology of the magnetic surfaces 
along which the wind flows. 

We also estimated the specific angular momentum extracted 
by the stellar winds from the star (lower panel in Fig. |4j defined 
as: 



Jsw 



Js 



w 



M sw 



(22) 



where the stellar wind torque 7sw has been obtained by inte- 
grating Eq. ( fTOt over the stellar surface from which the wind 
is launched, as in Eq. d21l . This last equation also provides 
the definition of the average magnetic lever arm Fa. The aver- 
age specific angular momentum of the case taken into account, 
jsw ~ 44 y/GM+R*, corresponds to a lever arm ?a ~ 21/?*. 
Notice that this value is slightly larger than the value found in 
case El (?a ~ 19/?*). Since the mass to magnetic flux ratio 
(77 = piip//?p) of the stellar winds of these two cases are com- 
parable, the different topology of the magnetic surfaces due to 
the interaction with the accretion funnels and the MEs likely de- 
termines the different lever arms. We just point out that the mag- 
netic configuration found in case El allows wide opening winds, 
while in case C03 the stellar wind assumes a more conical shape, 
where the presence of the MEs focuses the open magnetic flux 
towards the rotation axis. 

We can estimate the energy content of the stellar wind by 
evaluating Eq. (1201 at the stellar surface. It can be shown that 
this expression can be rewritten as: 



1- \2 



K* 



In 



y2 \rJ 2 vi 



(23) 



K* 



The first term in the right hand side, corresponding to the mag- 
netic flux, is the dominant one. Besides, the stellar boundary 
conditions on the outflowing material impose that h* * ^k*' 
so that the enthalpy provides the initial drive to the stellar wind: 
even if this term must be of different origin, it is just meant to 
mimic the effect of an extra pressure term. Therefore, we obtain 
that esw ~ 4.4V£ . This specific energy would correspond to 
a maximum asymptotic speed m PjCO as 2.9Vk*, but the outflow 
has attained a poloidal speed x y K * at the end of the computa- 
tional domain. The stellar wind has therefore the potential to be a 
light and very fast outflowing component, provided an efficient 
magnetic-to-kinetic energy conversion can be attained and this 
crucially depends on the asymptotic magnetic flux distribution. 
Notice that the total specific energy of the stellar wind is much 
higher than the one of the MEs (see Fig. [5j, almost one order of 
magnitude, mostly since the MEs are much more massive and 
therefore have less energy per particle available. This translates 
in a lower limit on the maximum terminal speed achievable by 
the MEs. 

We are not going to describe in great detail the properties 
of the disk wind accelerated along the open magnetic surface 



threading the accretion disk. Its mass outflow rate, obtained by 
integrating the mass flux equation Eq. ((9]) at the disk surface 
beyond R out , is rather small, see Fig. [3j and mostly concen- 
trated close to radius R out , where the magnetic field is stronger. 
Besides, it is clear from Fig. Q] that the Alfven surface lies very 
close to the disk: this means that the disk wind is characterized 
by a rather small Alfven radius, extracts a limited amount of an- 
gular momentum from the disk and therefore has a negligible 
impact on the angular momentum distribution of the star-disk 
system. This will be shown more quantitatively in the following 
Section. 



3.3. The disk angular momentum 

In this Section we analyze the impact of disk-driven outflows, 
MEs and disk winds, on the angular momentum distribution of 
the accretion disk. We start our analysis by taking into account 
the torques acting on the circumstellar disk and their effects on 
the disk dynamics. Accretion can be triggered (or hampered) by 
internal torques, e.g. turbulent viscosity, redistributing angular 
momentum radially inside the disk, or external torques, mag- 
netic stresses and mass transfers, extracting or supplying angu- 
lar momentum at the disk surface. This balance can be expressed 
in a steady situation by considering the angular momentum con- 
servation inside an annulus of the dis k of radi al width dr and 
thickness 2H. Following the notation of Paper j, the angular mo- 
mentum conservation can be formulated as: 



Face — Fi n t + F ma g + Tkin , 

where 

Face = — (^acc>" 2 £2disk) 



(24) 



(25) 



gives the angular momentum advection through the annulus. The 
mass accretion rate is defined as: 



Ma, 



pu r dz 

H 



(26) 



We express the internal torque F mt = F V1SC + Fb as the sum of the 
viscous "turbulent" torque: 



r ■ 

1 vise 



and the radial magnetic transport: 

„2 r+H 



_d_/V r +h 

dr\2j_ H 



B $ Brdz 



(27) 



(28) 



We included this last term for the sake of completeness: never- 
theless, the internal torque F; nt is dominated by the viscous term. 
The torque exerted by the large-scale magnetic field, extracting 
angular momentum at the disk surface, is defined as: 



r raag = sb;b; , 



(29) 



The kinetic torque, determined by the mass exchange at the disk 
surface, is given by: 



F km = -47Tr 3 p + Q. + U ; , 



(30) 



According to these definitions, a positive torque in the right hand 
side of Eq. ( 1241 extracts angular momentum from the annulus, 
thus favoring accretion. 
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Fig. 7. Upper panel. Logarithmic density map with sample field 
lines (yellow solid lines) and speed vectors (blue arrows) super- 
imposed. Middle panel. Radial profiles at the disk midplane of 
the rotation speed of the accretion disk Q/£2* {solid line), ac- 
cretion sonic Mach number {dot-dashed line) and plasma B {dot- 
ted line). The Keplerian and the Q = £2* rotation profiles are 
plotted with a dashed line. Lower panel. Radial profiles of the 
specific torques acting on the disk (see the text for definitions): 
magnetic (F mag , solid line), internal (r"j nt , dashed line and kinetic 
(r\i n , dot-dashed line) torques. In the three panels we marked 
with vertical lines the position of the truncation radius {R t ) and 
the anchoring radii at the disk surface of the outermost magnetic 
surface steadily connecting the star and the disk {R m ) and of the 
innermost open field line threading the disk {R ud- The panels 
represent temporal averages over 54 stellar rotation periods. 



In Fig.|7]we show the the right hand side torques of Eq. ( l24l i 
as a function of the radius r (lower panel), in the middle panel 
the effect of these torques on the disk structure (disk angular 
speed, accretion sonic Mach number, M s = \u r \ / -\/P/p\^ and 
plasma beta, B - 87rB 2 /P\ ), while in the upper panel we dis- 



play the corresponding density maps with field and stream lines 
superposed. The three panels are not snapshots at a given time, 
but have been obtained by averaging over 54 rotation periods 
of the star, from time t = 38 up to t — 92, in order to smooth 
out possible transient features. Besides, we marked three radii 
corresponding to the truncation radius {R t ), the anchoring ra- 
dius onto the disk surface of the outermost magnetic field line 
steadily connecting the disk with the star (/?;„) and of the inner- 
most open field line threading the disk {R ut)- These radii allow 
to distinguish three different zones: in the region between R t and 
Ri n , shaded in red in the middle panel, the disk can directly ex- 
change angular momentum with the star and form the accretion 
columns; the region between R m and R oa , shaded in green in the 
middle panels, despite being magnetically linked to the star, does 
not exchange angular momentum directly with it, since the mag- 
netic surfaces have expanded too much to be causally connected 
directly with the star. The material ejected from this region es- 
capes the stellar potential well instead of being accreted and the 
angular momentum extracted at the disk surface is transferred to 
the outflowing material: as already pointed out, this is the region 
from which the mass coming from the disk is accelerated to fuel 
the MEs. The disk region outside R out is threaded by open field 
lines along which, depending on the magnetic field strength, a 
disk wind can be accelerated. 

The lower panel shows that the internal turbulent transport is 
responsible for driving accretion on a large scale (r > 8 /?*). 
Even if the disk is threaded by the magnetic flux left by the 
opening of the dipolar magnetosphere, at this distance from the 
star the large-scale field is too weak to accelerate a powerful 
enough disk wind to exert a noticeable torque. Getting closer to 
the star, the large-scale open magnetic field threading the disk 
starts to be strong enough to influence the accretion dynam- 
ics: starting already in the region outside R oul , the magnetically- 
driven disk wind increasingly contributes to drive the accretion 
flow. Correspondingly, the disk dynamics start to change behav- 
ior (middle panels): the rotation profile is still Keplerian but, due 
to the growing magnetic torque, the accretion Mach number in- 
creases towards trans-sonic values. Since these outflows do not 
affect the disk angular momentum distribution in a relevant way 
(the rotation stays Keplerian) and extract a negligible amount of 
accreted mass, they will not be discussed in greater detail. 

The magnetic torque progressively becomes dominant in the 
disk region connected to the star from which the MEs arise 
{Ri n < R < R ont ). The middle panel shows that, in this region, the 
accretion Mach number starts to decrease after it has attained a 
maximum, almost sonic value: this corresponds to an adiabatic 
compression due to the push of the accretion flow against the 
stellar magnetosphere which is acting as a magnetic wall. This 
compression determines the enhanced mass-loading of the mag- 
netospheric ejections, as the kinetic tor que curv es show clearly 
in the lower panel. As already noticed in Paper I, this same effect 
is crucial to mass-load the accretion funnels. Besides, the ejec- 
tion torque (kinetic plus magnetic) becomes strong enough so 
that in this region the disk rotation becomes sub-Keplerian and 
even sub-stellar. In order to provide a more precise idea of the 
amount of angular momentum extracted by the MEs from the 
star-disk system, we can consider the angular momentum flux 
carried by the accretion flow through a vertical section of the 
disk at R out : 



J at 



= -2nr 



pU r U,p — T r( p — 



B,pB r 



dz 



(31) 



Without any other interaction of the disk with the surroundings 
inside R out , this would be the spin-up torque exerted by the ac- 
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cretion flow onto the star. We plot in the upper panel of Fig. 
3] the time evolution of the specific angular momentum carried 
through R om : 



Jacc,R m 



J-dcc,R.„ 



M, 



(32) 



acc,/? tll 



where the disk accretion rate M a cc,R rall has been obtained by eval- 
uating Eq. (l26i i at R out . The specific angular momentum y a cc,fi ollt 
is equal, on average, to y'acc,s OT , = 1-8 y[GM~*R~l. Notice that, 
since the definition of / a cc,s 0M (Eq. [3TJ includes the viscous and 
magnetic torques, j a cc,R om (Eq. 1321 is lower than the specific 
angular momentum expected from the advection process only 
(~ y/GM+R out m 2.6 y/GM+R*). Using the estimates done in 
Section 13.11 we can see that the MEs extract directly from the 
disk a fraction (M M E,d./ME,d) / (M acc , Rom ; acCrRoul ) * 53% of the an- 
gular momentum carried by the disk through R out , forcing the 
accretion disk to rotate at a sub-stellar rate. These ejections have 
therefore the important effect of extracting a relevant fraction 
of the disk angular momentum and consequently reducing the 
accretion torque. Besides, the material ejected from the disk to 
fuel the MEs clearly rotates slower than the star: since the mass 
ejected from the disk is also magnetically connected to the star 
and largely dominates the MEs inertia, it can extract angular mo- 
mentum from the star thanks to a differential rotation effect. 

The disk region inside R[ n can exchange its angular momen- 
tum directly with the star. Therefore, all the angular momentum 
that flows in across R[ n is eventually accreted by the star and thus 
determines a spin-up torque. In this magnetically connected re- 
gion, the stellar rotation tries to force the disk to co-rotate with 
it. Since in this example the disk rotates slower than the star at 
Ri n , the stellar rotation spins-up the disk back to £2*, thus exert- 
ing a negative magnetic torque on the disk (see the F mag curve 
in the lower panel of Fig. |7]i and extracting a fraction of the stel- 
lar angular momentum back to the disk. Notice that the r mag 
curve becomes positive again close to the truncation region and 
the disk starts to transfer angular momentum to the star. The 
truncation region is dominated by the kinetic torque, due to the 
mass loaded onto the base of the accretion funnelfl The exter- 
nal torque (magnetic plus kinetic) is characterized by a double- 
peaked profile, the inner positive peak being associated with the 
star-disk angular momentum exchange and the outer one with 
the torque exerted on the disk by the magnetospheric ejections. 

3.4. The stellar angular momentum 

In this Section we evaluate the impact of accretion and stellar 
outflows (MEs and stellar winds) on the temporal evolution of 
the angular momentum of the central star. We can express the 
time derivative of the stellar angular momentum /* by integrat- 
ing the angular momentum conservation equation over the stellar 
volume: 



A/* 
df 



-/ a cc,s + JmE,s + Jsw 



(33) 



where we separated the contributions to the torque due to ac- 
cretion, magnetospheric ejections and stellar winds. A positive 
(negative) angular momentum flux in the right hand side of 
Eq. (l33T > correponds to a spin-up (spin-down) torque. The three 



3 For a discussion on the dynamics of the truncation region and 
the acc retion funnels we refer the reader to Bessolaz et al. ( 2008|) and 
iPaperll (Section 3.1). 




Fig. 8. Temporal evolution of the torques acting directly onto the 
star, normalized to the stellar angular momentum. Plotted are 
the accretion torque (solid line), the stellar wind torque (dotted 
line) and the torque exerted by the MEs onto the stellar surface 
(dashed line). Conventionally a positive (negative) torque spins 
up (down) the stellar rotation. 



torques have been obtained by integrating Eq. (TTOb over three dif- 
ferent parts of the stellar surface: the area threaded by open field 
lines to evaluate the stellar wind torque 7sw (see Section [3~2l i: the 
area threaded by magnetic surfaces that undergo periodic infla- 
tion/reconnection events to evaluate the MEs contribution 7me,s 
(see Section [3~TT ); the area threaded by field lines steadily con- 
necting the star with the disk t o evalua ted the accretion torque 
Jacc.s- As already pointed out in iPaper II the contribution of the 
kinetic terms is completely negligible at the stellar surface, both 
for accreting and outflowing components, which are completely 
dominated by the magnetic torque. In the following discussion 
we will consider the total torques, keeping in mind that the mag- 
netic contribution is prevailing. In Fig. [8] we show the temporal 
evolution of the torques acting on the star, normalized to the stel- 
lar angular momentum = 7*Q*, where /* = tfM+R^ is the 
stellar moment of inertia. We assumed the typical normalized 
gyration radius of a fully convective star, i.e. k 2 = 0.2. Using 
this normalization, the curves in Fig. [8] provide directly the in- 
verse of the characteristic braking (or speed-up) timescale. To 
retrieve the physical units, the plots must be multiplied by Eq. 
(O that, expressed in terms of the stellar magnetic field inten- 
sity, takes the form: 



J_ 
J* 



= 10" 



(— ) 

UkG/ 



_M*_ 

0.5 Mr. 



-3/2 



R*_ 

2R C 



5/2 



yr 



(34) 



3.4.1. Accretion torque 

A rotating accretion disk is likely to exchange angular momen- 
tum besides mass onto the central object, thus providing a spin- 
up torque to the latter. The accretion spin-up torque can be 
parametrized as 7 acCjS = M acc s y acc s , where y acCjS is the specific 
angular momentum transported by the accretion streams and 
also along the magnetic surfaces connected to the disk which 
are not mass-loaded. A common parametrization for the specific 
accreted angular momentum is y acc = -\ZGM+R t , implying that 
a Keplerian accretion disk transfers to the star the angular mo- 
mentum possessed in the truncation region. In our simulations 
we can estimate the accuracy of this approximation. For exam- 
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pie the specific angul ar mom entum transferred by the disk to the 
star in case El from Paper I is larger than this reference value, 
jacc.s ~ 1.2 ylGM+Ri. since in this case the stellar magnetosphere 
is connected to the accretion disk over a large extent, even be- 
yond the corotation radius, the star can extract angular momen- 
tum from the disk in the entire region from R t up to R co . On 
the other hand, we already pointed out that in the "compact" 
magnetic configuration depicted in Fig. [1] the accretion torque 
is given by the angular momentum flux through a disk surface 
from R t up to R m < R co . In case C03, the average accreted an- 
gular momentum is approximately equal to 60% of the reference 
Keplerian value ylGM+R t . As already pointed out in Section [3~3l 
this effect is due the presence of MEs, which are extracting a 
consistent fraction of the angular momentum of the accretion 
disk before it is accreted onto the star. 



3.4.2. Stellar wind torque 

In Section [3721 we estimated the mass outflow rate and the spe- 
cific angular momentum extracted by the stellar wind from the 
star in the reference case C03. We can now compare the stel- 
lar wind torque with the accretion torque. On average, the stel- 
lar wind spin-down torque extracts around 66% of the accretion 
torque. This efficiency seems to be rather high, given the low 
ejection efficiency of the win d. For e xample, we recall that the 
stellar wind of case El from Paper I, was able to balance 20% 
of the accretion torque only. Three effects have enhanced the ef- 
ficiency of the spin-down torque: the mass ejection efficiency 
in case C03 is slightly higher than in case El (1.6% vs. 1.2%); 
in case C03 the magnetic lever arm is larger (?a ~ 217?* vs. 
7a ~ 19/?*). But the most important and interesting effect is due 
to the fact that in case C03 the MEs have already extracted a 
substantial amount of the disk angular momentum, reducing the 
accretion torque by a factor around 50% (see Section 13.31 1 and 
therefore enhancing the efficiency of the stellar wind torque. 



3.4.3. Magnetospheric ejections torque 

As we showed in Section 13.11 besides reducing the accretion 
spin-up torque, MEs are able to exchange angular momentum 
directly with the star. As we have already shown, the angular 
momentum exchange with the star is mainly controlled by the 
differential rotation between the star and the MEs: if at the cusp 
of the field line the mass loaded from the disk rotates slower 
(faster) than the star, the MEs exert a spin-down (spin-up) torque. 
Consistently, since we find that in case C03 MEs rotate slower 
than the star (see Section [3~3l l. they exert a net spin-down torque, 
equal to * 74% of the accretion torque. 

Summarizing, in the case C03 the combined action of MEs 
and stellar winds is able to balance the spin-up due to the ac- 
cretion torque. Notice that, even if MEs represent the dominant 
effect, extracting globally around 88% of the disk angular mo- 
mentum accreted through R out , the stellar wind can contribute 
significantly, extracting, in our example, around 3 1 % of the disk 
angular momentum flux through R out . Therefore, in the specific 
example analyzed in the previous Sections, the star is subject to 
a net spin-down torque: using Eq. ( 1341 to normalize to physi- 
cal units the sum of the three torques plotted in Fig. [8] we can 
estimate a characteristic spin-down timescale around 5.5 x 10 6 
yr. 
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Fig. 9. Same as Fig. [3] bur for case CI 



4. Varying the mass accretion rate 

In this Section we analyze cases CI and C01 and compare them 
with the results obtained from case C03, extensively presented 
in Section [3] We recall that these two cases are characterized by 
a different disk viscosity coefficient a v (see Table[TJ, with all the 
other parameters of the problem left unchanged. Primarily, we 
use the a v parameter to change the mass accretion rate of the 
disk. On the other hand, this modifies also the magnetic Prandtl 
number !P m = v v / v m of the disk. Therefore, bear in mind that our 
cases are not just characterized by different accretion rates and 
that other important aspects of the disk physics change. For ex- 
ample, the viscous accretion timescale varies from case to case, 
while the different Prandtl number can have an impact on the 
inclination of the field lines at the disk surface and on the ad- 
vection of the magnetic flux in the part of the disk dominated 
by the viscous torque. Different accretion rates could also have 
been obtained by varying the disk density with a fixed Prandtl 
number, possibly giving somewhat different results. 

4.1. High accretion rate and stellar spin-up 

Case CI is characterized by a value a v = 1, determining an ac- 
cretion rate around three times higher than the reference case 
C03, as confirmed by the temporal evolution of the accretion 
rate measured onto the stellar surface plotted in Fig. [9] There 
is another feature that is worth noticing: while in case C03 the 
accretion rate attains a rather constant asymptotic value, in case 
CI it slowly grows in time. We have already shown in Section 
13.31 that the accretion flow is controlled by the internal viscous 
torque far from the star, while it is determined by the external 
magnetic torques in the region of magnetospheric interaction. 
A mismatch between the two torques can trigger a long-term 
variability. For example, it seems that in case CI the magneto- 
sphere cannot sustain the accretion rate of the outer disk. This 
determines a density and pressure buildup inside the disk that 
causes, due to the a parametrization, an increase of the viscous 
torque and therefore of the accretion rate. The constant asymp- 
totic value of the accretion rate in case C03 is likely to be due to 
a better matching between the inner magnetic and outer viscous 
torques. 

The global magnetic configuration of case CI is very sim- 
ilar to case C03 (see Fig. Q]) and the same dynamical features 
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Fig. 10. Same as Fig. [7] but for case CI 



can be found. Regarding the MEs, their total mass outflow rate 
corresponds to around 6% of the disk accretion rate and, as in 
case C03, the mass extracted from the disk largely dominates the 
mass content of the MEs (see Fig . [9j . In the bottom panel of Fig. 
[TUlwe can see that the MEs are accelerated from a region of the 
disk located inside the Keplerian corotation radius which is more 
compact than in case C03, probably because of the larger mag- 
netic Prandtl number. Outside R oul , accretion is mainly triggered 
by the viscous torque, with a relatively small contribution from 
a weak disk wind, which, nevertheless, accelerates the accretion 
flow towards an almost sonic accretion speed (middle panel in 
Fig. [Tol l, without modifying the Keplerian structure of the disk. 
Inside the launching region of the MEs (R m < R < R out ), the 
torque due to the MEs becomes dominant. In order to estimate 
the efficiency of this torque we calculated the specific angu- 
lar momentum carried by the accretion flow through R out (see 
Eq. (|3TT >) and the specific angular momentum extracted by the 
MEs from the disk (see Eq. fFft ) and plotted them in the up- 
per and third panel in Fig. [TJJ respectively: MEs extract about 




Fig. 11. Same as Fig[4] but for case CI 



(M ME ,diME,d) / (M acc>Sou j acc , Roiil ) * 20% of the angular momen- 
tum carried by the disk through R oat , The lower torque efficiency 
measured in case CI is due to a lower mass ejection efficiency 
but also to a smaller specific angular momentum extracted from 
the disk: this is consistent with the fact that in case CI the mag- 
netic surfaces along which the MEs are accelerated have a larger 
inclination angle with respect to the disk surface (compare the 
upper panels of Fig.l7land[T0b. thus determining a smaller mag- 
netic lever arm. Nevertheless, the MEs torque is strong enough 
so that the disk rotation speed becomes sub-Keplerian (middle 
panel in Fig.fTOl, but not sub-stellar as in case C03. 

The angular momentum carried by the disk through is fi- 
nally transferred to the star and exerts a spin-up torque. In the 
lower panel of Fig. [10] it is possible to see that in the disk re- 
gion directly connected to the star below R[ n the angular mo- 
mentum is extracted from the disk and transferred to the star by a 
combination of magnetic and kinetic torques, where the peak of 
the latter corresponds to the angular momentum loaded onto the 
base of the accretion funnels. The external torque (magnetic plus 
kinetic) in the region R < R out presents qualitatively the same 
double-peaked profile of case C03, with the inner peak corre- 
sponding to the angular momentum exchange with the star, while 
the outer one representing the torque exerted by the MEs onto the 
disk. Using the same parametrization employed in Section [3.4.1l 
we can estimate that the specific angular momentum transferred 
to the star corresponds to /' acc s x 0.8 y[GM~*R t of the reference 
"Keplerian" value: this clearly shows that also in case CI one 
important effect of the MEs is to reduce the amount of angu- 
lar momentum that the disk transfers to the star, although less 
efficiently than in case C03. 
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Fig. 12. Same as Fig. [8] but for case CI 



The temporal evolution of the accretion torque is plotted 
in Fig. Q~2] together with the torques exerted by the MEs and 
the stellar wind onto the stellar surface. The asymptotic in- 
crease of the accretion torque is clearly linked to the temporal 
growth of the mass accretion rate discussed previously. The spin- 
down torque exerted by the stellar wind corresponds to 7.5% 
of the accretion torque only. Despite carrying a specific angu- 
lar momentum comparable and even slightly larger than case 
C03 (see bottom panels of Fig. [4l and fTTT>. the low ejection rate 
(Msw ~ 2 x 10~ 3 M acc , s ) strongly limits the efficiency of the spin- 
down torque. By inspecting Table |2l it is possible to notice that 
in case CI the connected magnetosphere contains more magnetic 
flux (<Jcd + ^mc) and it is even more compressed (see R cm ) than 
in case C03: the accretion funnels are located at slightly higher 
latitudes, leaving less room for the acceleration area of MEs and 
stellar winds. 

Besides, the stellar wind torque slowly drops in time due 
to the decrease of the stellar wind specific angular momentum 
(lower panel in Fig. fTTb : this is likely associated with the in- 
crease of the mass accretion rate which can modify the shape of 
the magnetic nozzle at the base of the stellar wind, the acceler- 
ation efficiency of the outflow and therefore its magnetic lever 
arm. The torque exerted by the MEs directly onto the star is al- 
most completely negligible. Besides, it is possible to see in Fig. 
[12] that for t > 60 it even becomes positive, contributing to the 
spin-up torque, even if to a very small extent. As we pointed 
out in Section [3.4.3l the angular momentum exchange between 
the MEs and the star depends on a differential rotation effect be- 
tween the star and the mass loaded from the disk onto the MEs. 
We saw in Fig.[10]that in case CI the base of the MEs rotate at a 
sub-Keplerian but still super-stellar angular speed: consistently, 
the MEs do not directly brake the stellar rotation and can actu- 
ally spin it up. The same effect is also shown in the two central 
panels of Fig. QT| the total specific angular momentum carried 
by the MEs (second panel) is smaller than the specific angular 
momentum extracted from the disk (third panel). In case CI, 
contrary to case C03, the MEs can extract energy and angular 
momentum from the disk but then transfer a small fraction back 
to the star. 

Since in this case the spin-down torques are very inefficient, 
the star experiences a strong spin-up torque due to accretion. 
Normalizing the torques plotted in Fig. Q~2] with Eq. ( 1341 1. we 
estimate a characteristic spin-up timescale around 5.5 x 10 5 yr. 



4.2. Low accretion rate: 

transition to a "propeller" regime 

Case C01 is characterized by a lower disk turbulent viscos- 
ity (a v = 0.1) and therefore a lower mass accretion rate. 
Consistently with the results of cases C03 and CI, in case C01 
the viscous torque likely supports an accretion rate which is 
smaller than the one determined by the inner magnetospheric 
torques. Therefore, the disk tends to empty, reducing its surface 
density and pressure, thus decreasing the a torque and its accre- 
tion rate. Besides, the a-disk model that we employed as initial 
condition presents a large-scale meridional circulation pattern 
for low ff v values, with the disk accreting along the surface lay- 
ers only and excreting along the midplane. Specially for low a v 
values, this model becomes very sensitive to local changes of 
thermal pressure gradients and its accretion rate is difficult to 
control. 

Due to these effects, the disk accretion rate in this case is 
rather low and slowly decreases in time. This is clearly visible 
in the left panel of Fig. Q~3] where, after a strong initial tran- 
sient peak, the accretion rate diminishes steadily, subsequently 
it starts to oscillate and, at last (f > 60), it becomes highly un- 
steady and intermittent, varying periodically between relatively 
high and extremely low values. The system evolution during one 
of these cycles is depicted in Fig. [14] which clearly shows how 
the accretion cycles correspond to a periodic oscillation of the 
truncation radius. 

The low-accretion phases correspond to the truncation radius 
moving close to the Keplerian corotation radius: in this situation 
it is hard to form the accretion funnels, since the centrifugal bar- 
rier raised by the rotating magnetosphere prevents the disk ma- 
terial from falling towards the star and an extra th ermal pressure 
gradi ent would be required to cross the barrier (|Koldoba et al.l 
2002). As extensively discussed by iBessolaz et alJ d2008h . the 
disk truncation is primarily determined by a pressure equilib- 
rium between the disk and the magnetosphere. Therefore in case 
C01, as the accretion rate decreases, the disk reduces its push 
against the magnetosphere and the truncation radius progres- 
sively moves outwards closer to corotation. A consequence of 
the disk being truncated beyond the corotation radius is that, 
since in this region a Keplerian disk rotates slower than the 
star, the stellar rotation tries to increase ("propel") the disk an- 
gular speed in the region directly connected to the star, hence 
the appellative "prope ller" regime (llllarionov & Sunvaevll 19751 
Ustvugova et al] |2006f) . Notice that, in our simulations, a pro- 
peller effect can be present even below the Keplerian corota- 
tion radius, whenever the disk rotation becomes sub-stellar, see 
for example Fig. [7] Despite the presence of a propeller effect, 
the angular momentum extracted from the star is likely to be 
energet ically insufficient to g ravitationally unbind all accreting 
matter ( I Spruit & Taaml [19931). Besides, as the disk moves far- 
ther from the star, the stellar magnetic field weakens and the 
magnetosphere gets more and more inflated, losing part of its 
connection with the disk. As the field lines open up, a massive 
magnetospheric outflow is launched along the opening magnetic 
surfaces. This ejections are analogous to the MEs described in 
the previous Sections; in a propeller phase, they can be even 
stronger and can extract an important fraction of the angular 
momentum of the disk. Because of the torque exerted by these 
strong ejections, the disk temporarily increases its accretion rate 
and its drive against the stellar magnetosphere. The truncation 
radius can therefore move closer to the star, the funnel flows are 
fueled again and the accretion rate onto the surface of the star 
increases. At the same time, the inflated field lines along which 
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Fig. 13. Left panel. Temporal evolution of different mass fluxes in case C01. Plotted are the mass accretion rate onto the stellar 
surface {solid line), the stellar wind outflow rate (dashed grey line) and the stellar mass flux contribution to the MEs (dashed line). 
Right panel. Temporal evolution of the torques acting onto the stellar surface, normalized to the stellar angular momentum. Plotted 
are the spin-up torque associated with accretion (solid line), the spin-down torque exerted along the magnetic surfaces connecting 
the star and the disk (dashed line) and the spin-down torque due to the stellar wind (dotted line). 



the MEs have been launched tend to close again and quench the 
outflow (see the lower central panel in Fig.[T4b. Instead of losing 
angular momentum to the outflow, the disk can now acquire it 
from the star along the field lines that are connected to the star, 
since, due to the action of the MEs, the disk rotation is at least 
in part sub-stellar, analogously to case C03. The disk therefore 
decreases its drive against the compressed magnetosphere which 
pushes back the disk towards the corotation region and the cycle 
repeats. Notice that the disk can move outwards not only if it be- 
comes super-Keplerian, but, since during the accretion phase the 
closed magnetosphere is strongly compressed (see for example 
Fig. 8 in lPaper m . the poloidal magnetic pressure can effectively 
push the disk outwards whenever the latter reduces its thrust. 

These accretion/ejec tion cycles have been observed in other 
nume rical works (e.g. iGoodson et al.l Il999bt lUstvugov a et al.l 
2006): notice that the range of timescales and amplitudes of the 
oscillations observed in these papers can be simply due to dif- 
ferent parameters of the models, i.e. stellar rotation period, mag- 
netic field intensity, disk accretion rate. The accretion cycles of 
a disk truncated clos e to corotati on have been predicted a lso by 
ISpruit& Taaml(ll993h (see also lD'Angelo & Spruitj|2010l) . They 
showed how an accretion disk truncated beyond the Keplerian 
corotation radius can possibly readjust its temperature and den- 
sity structure so that the viscous stresses get rid of the excess 
angular momentum coming from the star and the accretion flow 
can cross the corotation radius and cyclically form the accretion 
funnels. The characteristic period of their cycles is obviously as- 
sociated with the viscous accretion timescale in the corotation 
region. In our simulation the timescale of the accretion cycles 
is much shorter: since the oscillations are mainly driven by the 
torque exerted by the MEs, the typical period is of the order of a 
few Keplerian orbits around the corotation radius. 

The torques exerted onto the star are plotted in the right panel 
of Fig. Q~3] Clearly, during the "propeller" phases characterized 
by an extremely low accretion rate, there is no spin-up torque 
associated with accretion. Besides, a strong spin-down torque is 
exerted along the field lines connected to the disk. Because of 
the large amount of angular momentum extracted by the magne- 
tospheric ejections, both the accretion disk and the MEs rotate 



slower than the central star, even in the sub-corotation region. 
Therefore the star can be efficiently spun-down along the mag- 
netic surfaces directly connected with the disk and the MEs. 

During the high accretion phases, the disk can deposit its an- 
gular momentum along the funnel flows, exerting a small but 
noticeable spin-up torque. The spin-down torque exerted along 
the closed magnetosphere is reduced, since the MEs, which are 
the main responsible for this torque, are weaker during the accre- 
tion phases: nevertheless, it is still possible to balance the torque 
due to accretion. The accreting phases of case C01 resemble, at 
least qualitatively, the steadily accreting case C03. 

Finally, it is important to notice that also in case C01 a strong 
stellar wind is present, exerting an important spin-down torque 
onto the star. Table|2]shows that in this case the stellar wind can 
exploit a larger amount of open stellar flux that becomes even 
more important during the non-accreting phases, when the stellar 
wind torque seems to become even stronger. This is obviously 
not consistent with having a stellar wind fueled by the accretion 
power: in case C01 the stellar wind would require a considerable 
driving power, as well as during the phases during which the disk 
is not accreting. On the other hand, in this case the role played 
by the stellar wind can be neglected, since the torque exerted 
by the star-disk-MEs interaction is largely sufficient to brake the 
stellar rotation. Even neglecting the stellar wind torque, we can 
estimate that, on average, the characteristic spin-down timescale 
in case C01 is approximately equal to 8 x 10 5 yr. 

5. Discussion 

In this Section we compare our findings with works that have 
described related scenarios and addressed similar issues. First of 
all, we point out that phenomena qualitatively similar to the MEs 
have been observed in a number of simulations of magnetic star- 
disk interaction: in spite of differences in the interpretation of 
the results, we have the feeling that these numeric al experiments 
show essentially the same process. For example. lHavashi et al.l 
(119961) observed a plasmoid ejection associated with the ex- 
pansion of the magnetosphere and associated the reconnection 
episodes to the X-ray flares observed in young stars. A series 
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Fig. 14. Time evolution of an accretion cycle during the propeller phase of case C01. We show logarithmic density maps with 
sample field lines (solid lines) and speed vectors (blue arrows) superimposed. The yellow solid line follows the evolution of a single 
magnetic surface, showing clearly the periodicity of the accretion-ejection cycles. Time is given in units of the stellar rotation period. 



of papers dGoodson et al.ll 19971 Il999allbt iMatt et alj|200l pro- 
posed the "epis odic magnetospheric inflation" mechanism (EMI, 
IMatt et ai1l2003l) to explain the origin of jets from young stars: 
this is probably due to the fact that, besides the presence of 
uncollimated outflows commonly related to this phenomenon, 
these simulations show the formation of a dense axial collimated 
jet, fueled during the inflation phases by mass coming from the 
accretion streams and focused around the rotation axis by the 
expansion of the closed magnetosphere. While their description 
of this mechanism is very phenomenological, this dense col- 
limated axial feature has not been observed in our numerical 
models. Additionally, the periodic behavior of the mass accre- 
tion rate suggests that t hese solutions m i ght be in a "propeller" 
regime. More recently, iRomanova et al.l (120091) developed both 
axisymmetric and three dimensional non-axisymmetric simula- 
tions that describe the formation of uncollimated outflows ("con- 
ical winds") coming from the boundary of the closed magneto- 
sphere, both for slowly and fast rotating stars. Despite clear sim- 
ilarities between these results and the scenario proposed in this 
paper, we have performed a much more detailed analysis of the 
dynamics and the energetics of MEs, which allowed us to obtain 
precise results about how they can affect the angular momen- 
tum of the star-disk system. For example, we showed how MEs 
can extract angular momentum both from the disk and the star 
so that their action can be compared to that of X- Winds on one 
hand and stellar winds on the other. 

Similarly to X- Winds, we showed that MEs can extract a rel- 
evant fraction of the disk angular momentum in order to reduce 
the spin-up torque exerted by accretion. Actually, the X-Wind 
model represents an extreme solution in which all the angular 
momentum carried by the accretion flow is extracted from the 
disk, so as to cancel at least the accretion torque. For the MEs 
we estimated a lower efficiency, e.g. 53% in the fiducial C03 
case. Due to the lower torque efficiency, the energy needed to 
drive MEs is compatible with the mechanical power released by 
the accretion flow, while t he latter is likely to be i nsufficient to 
drive an efficient X-Wind dFerreira & Cassdl2012l) . In addition, 
the MEs can also extract angular momentum directly from the 
star, similarly to stellar winds. The main difference is that, while 
the spin-down torque exerted by stellar winds strong ly depends 
on the mass flux of the wind (Matt & Pudritz 2008a), the torque 
of the MEs is due to a differential rotation effect between the star 



and the material ejected from the disk by MEs. Therefore, MEs 
do not seem to suffer the energetic problems associated with the 
mass-loading of the stellar winds. Notice that the idea of having 
an outflow mass-loaded from the accretion disk capable of ex- 
tracting angular momentum fr om the central star i s close to the 
Reconnection X-Wind model dFerreira et al.ll2Q00b . even if this 
scenario had been envisaged for a different magnetic configura- 
tion. 

An important difference between MEs and stationary winds 
is represented by their asymptotic behavior. While the large- 
scale acceleration and collimation properties of steady outflows, 
whether coming from the star or the disk, depend on the global 
distribution of their poloidal magnetic flux and electric current, 
we showed that MEs rapidly disconnect from the central re- 
gion of the disk-star system due to reconnection events and 
propagate ballistically afterwards, without accelerating any fur- 
ther. Flowing in between a stellar wind and a disk wind (see 
Fig. [TJ, the confinement of the MEs depends on the collimat- 
ing/decollimating behavior of the other two outflows (see dis- 
cussion in AppendixlAl. 

A final comparison can be made between the poloidal mag- 
netic configuration found in our simulations and the one ex- 
pected from th e X-Wind model, at le ast in its "dipolar" formu- 
lation (see e.g. lOstriker & Shul[l995l) . The amount of magnetic 
flux contained in the magnetic cavity is rather consistent with 
our results (<£>mc in Table The X-Wi nd model also predicts 
that the remaining flux (see Eq. (11) in Mohan tv~& Shul f2008) 
is trapped in the disk in a small region around the corotation ra- 
dius and is evenly shared among the accretion funnels (4>acX 
the open flux (<t>sw) and the remaining closed field lines which 
do not accrete (Ome + ^cd - 3*ac)- In our simulations the flux 
repartition can noticeably vary from case to case but it can be 
actually similar to the X-Wind distribution in cases close to the 
spin equilibrium (e.g. during the accretion phases of case C01). 
Three important differences can be pointed out: (1) in our simu- 
lations the magnetic flux involved in the star-disk interaction is 
distributed over a sizeable region of the disk, as shown by the 
anchoring radii of the different magnetic surfaces given in Table 
|2j (2) there is no magnetically "dead zone", namely a magnetic 
zone with neither mass nor angular momentum exchanges, as 
proposed in the X-wind scenario: in contrast, we obtain an ex- 
tremely active, outbursting zone where MEs take place; (3) the 



16 



C. Zanni and J. Ferreira: MHD simulations of accretion onto a dipolar magnetosphere 



steadily open magnetic surfaces threading the disk, where the 
X-Wind should be accelerated, play a marginal role in our solu- 
tions. 

In the following subsection we are going to compare more 
quantitatively the efficiency of the spin-down torques found in 
our simulations with the outcome of other scenarios proposed in 
the literature. 



5.1. Zero-torque condition 

By progressively lowering the disk accretion rate, the three sim- 
ulations analyzed in the paper show a transition from a strong 
spin-up (case CI) to an efficient spin-down state (case C01). 
This trend corresponds also to a different repartition of the stel- 
lar magnetic flux among the spin-down/spin-up phenomena (see 
Table 0, where the magnetic flux associated with the accretion 
spin-up torque (Oac and OcdX decreases at the expenses of the 
spin-down mechanisms (<1>me and <£>sw)- In case C03 the torques 
exerted by the MEs, with a weaker but non-negligible contribu- 
tion from a light stellar wind, are able to balance the accretion 
torque, so that the net total torque is approximately zero (actually 
slightly negative in this specific case). Correspondingly, the trun- 
cation radius moves from R t w 0.45/? CO (case CI) to R t w 0.8/? co 
(case C01), with the zero-torque configuration approximately lo- 
cated at R t aa Q.6R C0 . Besides, the sim ulations confirmed the re- 
sults of iBessolaz et al.l (|2008) and of Paperjl s howing that the 
Alfve n radius for a spherical free-fall collapse dElsner & Lambl 
Il977l) . namely 
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is a good parametrization of the truncation radius. In our sim- 
ulations we have found R t as OAR a- Therefore the zero-torque 
state corresponds roughly to a situation where R co as 0.67^a- A 
proportionality relationship between R co and R& has been used 
in different scenarios to estimate the stellar rotation period cor- 
responding to a zero-torque situation: 
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(36) 



where K = (Rco//?a) • Obviously, the longer the equilibrium 
rotation period, the more efficient the spin-down mechanism. We 
already showed that our simulations suggest a value K as 0.54. 
We can compare this result with the outcome of other popular 

scenarios. 

Applying the classica l jGhosh & Lambl scenario to the case 
of T Tauri stars, iKonigll dl99ll) used a value K = 0.87. On 
the other h and, the correct io ns brou ght to the iGhosh & Lambl 
model by iMatt & Pudritzl (l2005al) provide an upper limit 
K < 0.3, kee ping in mind that, according to iPaper J. the 
Ghosh & Lamb spi n-down torq ue is likely to be even weaker 
than the lMatt & Pudritzl d2005al) es timate. In the case of the X- 
Wind model, Ostrik er & Shul d 19951) employed a value K = 0.89, 
even if we recall that many properties of the scenario are based 
on ad-hoc assumptions and are not the result of a detailed dy- 
namical calculation. The results of the time-dependent axisym- 
metric simulations bv lLong et al.l (120051) have been recently re- 
examined, suggesting a value K = 0.52 (M. Romanova, priv. 
comm.), in good agreement with our findings. 



For the accretion powered stellar wind scenario it is possi- 
ble to derive an expression analogous to Eq. (l36*l l by equating 
the spin-down torque of a stellar wind (Jsw = A^sw^a^*) w i m 
the accretion torque of the di sk (7 acc oc M acc -\/GM^R^. Using 
Eq. (12) from lMatt& Pudritzl (|2008a) to determine the average 
Alfven radius ?a and assuming the accretion torque found in 
our case El since it is not affected by the presence of the MEs 
(y acc = 1.2M acc y/GM+R t with R t = 0.4fl A ), we obtain 
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(37) 



similarly to Eq. (17) in IMatt & Pudritzl d2008bi) . This equation 
clearly shows that an ejection efficiency M^/M- dcc ~ 0.1 is 
needed to give a characteristic rotation per i od com parable to the 
other scenarios, as discussed in lMatt et alj 

Since Eq. d36"b and d37l > describe an equilibrium between ac- 
cretion and spin-down torques, the same zero-torque condition 
(i.e. the same P eq ) can be obtained by varying accordingly the 
mass accretion rate, proportional to the accretion torque, and the 
stellar magnetic field intensity, mainly related to the spin-down 
torque. This implies that stars with a comparatively lower ac- 
cretion rate need a weaker dipolar field to balance the accretion 
torque. 



(35) 5.2. Stellar contraction and spin-down 



The zero-torque condition discussed in the previous Section is 
rather challenging to obtain for all the discussed models, but it 
is obviously not sufficient to keep constant the rotation period 
of a contracting star. The time derivative of the stellar angular 
velocity 
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depends also on the temporal evolution of the stellar mass and ra- 
dius. The right hand side of Eq. (l38l defines three timescales: the 
term M*/M* ~ l/f acc is associated with the accretion timescale, 
typically f acc ~ 10 7 - 10 10 yr for a CTTS; the term ~ 
-IAkh is associated with the the Kelvin-Helmholtz contraction 
timescale: 



?kh = 1.8 x 10 6 
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Clearly, it is necessary to exert a net negative torque onto the 
star with a characteristic spin-down timescale < at least 

comparable to the contraction timescale to maintain a steady 
stellar rotation period. In our simulations this condition is clearly 
satisfied during the propeller phases of case C01. This regime 
occurs whenever the disk is truncated close enough to the coro- 
tation radius (R t > 0.8/? co in our simulations). Notice however 
that this condition is necessary but not sufficient to have a spin- 
down torque strong enough to balance the stellar contraction. 
While different fi* -M acc combinations can satisfy the condition 
R ( oi 0.8/? co , as in the case of a null torque configuration, a strong 
magnetic field is necessary to provide a spin-down torque able 
to balance the stellar contraction. This is clearly shown by Eq. 
(1341 . used to normalize the torques plotted in Fig. 151 IT21 and IT31 
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In case C01 we must assume a strong dipolar component, typi- 
cally of the order of the kG, to obtain a short enough spin-down 
timescale. This result shows that, even if stars characterized by 
a low accretion rate (e.g. < 1O~ 9 M yr _1 ) need a dipolar field 
weaker than a kG to cancel the accretion torque, they still need 
a kG dipole to balance their contraction. 



6. Summary and conclusions 

In this paper we presented axisymmetric MHD time-dependent 
simulations of the interaction of a dipolar stellar magnetosphere 
with a surrounding viscous and resistive accretion disk (using a 
prescriptions). We assumed a magnetic coupling strong enough 
(i.e. a resistivity sufficiently low) so that the buildup of the 
toroidal field magnetic pressure due to the star-disk differen- 
tial rotation inflates and opens up the dipolar structure close 
to the central star and the truncation region. In particular, the 
strong coupling prevents the closed magnetosphe re to extend be- 
yond the Keplerian corotatio n radius , so that the lGhosh & Lambl 
spin-down model, studied in lPaper 1 cannot be applied anymore. 
On the other hand, the simulations showed that magnetospheric 
ejections naturally arise because of the process of expansion and 
reconnection of the magnetospheric field lines connected to the 
disk. We extensively studied the dynamical properties of these 
ejections with a special focus on their impact on the angular 
momentum of the star-disk system. At the same time we have 
included in our models the spin-down torque exerted by a stel- 
lar wind. We here summarize the main results of our numerical 
experiments: 

1. MEs can exchange mass, energy and angular momentum 
both with the star and the disk. Their inertia and mass content 
is largely dominated by the material loaded from the accre- 
tion disk. If the mass load rotates slower than the star the 
ejections can be powered by both the stellar and disk rota- 
tion, as in a huge magnetic slingshot. 

2. MEs cannot explain the jet phenomenon in T Tauri stars: 
(1) their terminal speed is unlikely to be higher than the 
gravitational escape speed; (2) after they disconnect from 
the star-disk system they propagate ballistically as uncolli- 
mated plasmoids. Their confinement depends on the colli- 
mation properties of the outflows between which they prop- 
agate, stellar and disk winds. 

3. MEs crucially contribute to control the stellar rotation pe- 
riod. On one hand, they efficiently extract angular momen- 
tum from the disk close to the truncation region so that the 
spin-up accretion torque is sensibly reduced. On the other 
hand, if the torque exerted onto the disk is strong enough so 
that the ejected plasma rotates slower than the star, MEs can 
extract angular momentum directly from the star thanks to a 
differential rotation effect (slingshot effect). We found a bal- 
ance between the accretion torque and the spin-down torque 
for R t * 0.6R co . 

4. The efficiency of the spin-down torque of MEs is comparable 
to other scenarios propos e d in th e literature, keeping in mind 
that: the lGhosh & Lambl d 19791) spin-down torque is highly 
overe stimated (see [Matt & Pudritz 2005al IZanni & Ferreiral 
2009); a fully self-consistent dynamical model of the X- 
Wind scenario is currently not available; the mass ejection 
efficiency of a stellar wind capable to balance at least the 
accretion torque is energetically very demanding. While the 
torques exerted by MEs onto the disk and the star share some 
similarities with the X-Wind and stellar wind models respec- 



tively, MEs do not seem to be affected by the energetic limi- 
tations that concern the other two scenarios. 

5. We limited the mass ejection efficien cy of the stellar winds to 
a few percent. Consistently with the Matt & Pudritzl (l2008bl) 
results, these winds are not sufficient to balance the spin-up 
due to accretion and contraction. Nevertheless, we found that 
for a wind mass flux around 1 - 2% of the accretion rate, the 
spin-down torque corresponds to 20 - 30% of the accretion 
torque. We found that a light disk wind is launched along 
the open magnetic surfaces threading the disk but, due to the 
weakness of the field, the impact of this outflow on the angu- 
lar momentum structure of the Keplerian disk is negligible. 

6. In a propeller phase, when the truncation radius gets suf- 
ficiently close to corotation (R t > 0.8/? co ), the spin-down 
torque exerted by the disk and the MEs can even balance 
the spin-up due to contraction. On the other hand, during the 
propeller phases the accretion becomes intermittent on a dy- 
namical timescale, corresponding to a few rotation periods 
of the star. Even if this effect is enhanced by the axial sym- 
metry of our models, there is no observational evidence of 
such a behavior. 

We want to conclude by pointing out a limit that affects all 
the scenarios discussed in this work, including the model pre- 
sented in this paper. All the scenarios are based on axisymmetric 
models of purely dipolar stellar magnetospheres and come to 
the conclusion that, for typical slowly rotating CTTS, a dipo- 
lar component around ~ 1 kG is needed to balance, at least, 
the spin-up torque due to accretion. Whil e this has been ob - 
served, for example, in the case of AA Tau dDonati et al.l20 10b). 
such a strong dipolar component d oes not seem to be a com- 
mon f eatu re among T Tauri s tars ((Donati et al] 120071 1201 Oal 
12011333: iHussainet al.l f2009). While a larger sample of mag- 
netic field measurements is clearly needed, stellar torque models 
must start to consider more realistic magnetic field configura- 
tions, in cluding non axisymmetric and multipolar components 
(see e .g. lLong et al.ll201l[ lRomanovaetalJl201lL IVidotto et all 
BOTlh . 
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Appendix A: Global magnetic coupling 

Since the star-disk interaction is mainly controlled by magnetic 
stresses, the best way to have a general view of the magnetic 
coupling between the different parts of the system, the accretion 
disk, the star and the outflows, is to look at the poloidal electric 
currents J p , circulating along the rB^ = const, isosurfaces. The 



poloidal component of the Lorentz force J p x B^ is perpendicu- 
lar to these surfaces, while the toroidal force J p xB p is obviously 
proportional to the magnetic torque. Besides, in the region close 
to the star, the dominant components of the angular momentum 
and energy poloidal fluxes are directed along the poloidal mag- 
netic field and are proportional to oc -B^,B p : the sign of gives 
therefore an information about the direction of the Poynting flux. 
In Fig. [15] we plot sample rB^ = const, isosurfaces to outline the 
poloidal current circuits characterizing cases CI and C03. The 
two panels of Fig. [16] refer to the current configuration during 
the high accretion (left panel) and low accretion phases (right 
panel) of the propeller stage (? > 60) of case C01. We can dis- 
tinguish three main current circuits labeled as A, B and C in the 
four panels. In circuits A and C the currents circulates clock- 
wise, while counterclockwise in circuit B. Different colors cor- 
respond to negative or positive values of Bf. dark-red circuits 
(B<p > 0) are responsible for extracting angular momentum and 
energy from the disk, while the light-yellow circuits (B^ < 0) 
extract angular momentum and energy from the star. 

The electromotive force of A-labeled circuits is due to the 
star-disk differential rotation: the current flows out from the 
stellar surface towards the disk, flows back along the accretion 
funnels and closes inside the star. This current circuit exerts a 
toroidal braking force inside the disk and the accretion columns 
and a spinning-up force inside the star, thus transferring angular 
momentum from the disk and the funnels to the star. It is there- 
fore responsible for the accretion spin-up torque plotted in Fig. 
[3]and[9] Notice how this circuit becomes smaller and smaller as 
the accretion torque decreases going from case CI to the high ac- 
cretion phases of case C01 and completely disappears in the low 
accretion stages of case C01, during which the accretion torque 
is completely negligible. 

Circuit C brakes down the disk rotation and is responsible 
for the magnetic driving of outflows launched from the disk, in 
particular for the magnetic acceleration of the part of the MEs 
mass-loaded from the disk. In the launching region of the MEs, 
it provides a strong vertical force which uplifts matter at the disk 
surface (see Fig. [6), thus contributing to the high mass ejection 
efficiency of these outflows. This circuit corresponds to the in- 
nermost part of the butterfly-s haped current ci rcuits characteris- 
tic of extended disk winds, see lFerreir al (fT997h . 

Current circuit B is associated with the energy and angular 
momentum extraction from the star. The current flows out from 
the stellar surface at mid latitudes and flows back to the star at 
higher latitudes. In Case CI the current flows mostly along the 
open field lines anchored onto the stellar surface, therefore fu- 
eling the stellar wind. There is no spin-down circuit associated 
with the MEs or the disk inside the closed magnetosphere. On 
the other hand, the other cases show clearly that the spin-down 
circuit B couples the star with the disk, the MEs and the stellar 
wind. These three dynamical constituents can in fact extract a 
fraction of the stellar angular momentum. The current flowing 
inside the disk, the MEs and the stellar wind provides in fact a 
J x B force which tries to spin-up the material rotating at sub- 
stellar rates, so that the star loses its angular momentum. The 
circuit B clearly identifies the parts of the disk still magnetically 
connected to the star which rotate slower than the latter. 

Altogether, the four panels show how the spin-up effects 
weaken and the spin-down action strengthens going from case 
CI to case C01. We can actually see that in case CI the MEs, 
which are efficiently extracting angular momentum from the 
disk, are transferring energy and angular momentum to the star, 
thus providing a spin-up torque, as we saw in Section l4~Tl In case 
C03 MEs extract energy and angular momentum both from the 
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Fig. 15. Poloidal current circuits flowing in the star-disk system in cases CI (left panel) and C03 (right panel). Dark (red) circuits 
circulate clockwise along isosurfaces of rB^ > 0. Yellow (light) circuits circulate counterclockwise along isosurfaces of rB<f, < 0. 
Current circuits are superimposed to a logarithmic density map. The figures have been obtained by time averaging over 54 stellar 
rotation periods. 




Fig. 16. Poloidal current circuits flowing in the star-disk system in case C01. Color codes are the same as in Fig. [15] The left panel 
refers to the accretion phases, while the right panel represents the strictly propeller phases. The left (right) panel has been obtained 
by averaging in time current and density during the maxima (minima) of the accretion rate for t > 60, see Fig.[T3l 



star and the disk and the two fluxes converge at the cusp of the 
field line. In the propeller phases of case C01, MEs seem to be 
powered almost exclusively by the stellar rotation. 

Finally, these figures clearly show that on a large scale MEs 
tend to propagate along the current sheet at the interface be- 
tween circuits B and C. While circuit B provides a decollimat- 
ing Lorentz force, circuit C tries to collimate towards the axis. 
Therefore, since MEs propagate ballistically after they have dis- 
connected, their collimation properties depend on the equilib- 
rium between the decollimating inner field lines and the colli- 
mating outer ones, i.e. on the equilibrium between the inner stel- 



lar wind and an outer disk wind. A proper treatment of this is- 
sue therefore requires simulations covering larger spatial scales, 
properly treating all the outflows present in the system. For ex- 
ample, the large-scale simulations presented by iGoodson et al.l 
dl999al) suggest that an inertial confinement due to a very thick 
disk and a rather dense corona can help focusing the closed mag- 
netosphere expansion and the propagati on of the magnetospheric 
outflows towards the axis (see also Li et al. 2001). Matt et al. 
d2003l) showed how a strong enough outer poloidal magnetic 
field anc h ored i n the disk can confine the EMIs. More recently, 
iLii et a D (120 12h showed that, at least for high mass accretion 
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rates, more typical of EXORs and FUORs, the innermost open 
disk field lines can be sufficiently mass loaded so that the corre- 
sponding disk wind can collimate the inner outflows. 



